DISPERSION RELATION FOR WATER WAVES 
WITH NON-CONSTANT VORTICITY 



PASCHALIS KARAGEORGIS 



Abstract. We derive the dispersion relation for linearized small- amplitude gravity waves 
for various choices of non-constant vorticity. To the best of our knowledge, this relation is 
only known explicitly in the case of constant vorticity. We provide a wide range of examples 
including polynomial, exponential, trigonometric and hyperbolic vorticity functions. 



1. Introduction 

Consider the classical water wave problem with a free surface over a flat bottom under 
the influence of gravity. Using bifurcation and degree theory, Constantin and Strauss pQ 
proved the existence of large-amplitude, two-dimensional, inviscid, periodic traveling waves 
with vorticity. A relation that arises in their work, also known as the dispersion relation, 
is the necessary and sufficient condition for local bifurcation. It is known explicitly in the 
case of constant vorticity [TJ H] but not for any other vorticity functions. Our goal in this 
paper is to provide a wide range of examples of non-constant vorticity functions for which 
the dispersion relation can be determined explicitly. 

Before we go on, however, let us first summarize the existence result of Constantin and 
Strauss PQ. For the sake of brevity, we only describe the details which are most relevant to 
our work and refer the reader to [T] for more details. We let (u, v) denote the velocity field 
and we assume that u < c, where c is the wave speed. This assumption makes the relative 
mass flux p negative, and it also allows us to express the vorticity u = j(if>) in terms of the 
stream function ■0. As a preliminary result to local bifurcation, we first recall the existence 
of a curve of trivial solutions (parallel shear flows), which are gravity waves with the bottom 
at y = and the top at y = d. Following [1], let us also introduce the notation 

r(p) = I l(s)ds, r min = minl» < 0. (1.1) 

JO \Po>°] 

For a proof of the following two results, we refer the reader to [TJ Section 3.1]. 
Lemma 1 (Parallel shear flows). For each A > — 2r m ; n , there exists a trivial solution 

r ds 

V(P) = / A+9TV \ (L2) 

which corresponds to a parallel shear flow with velocity field 

u(y) = c-y/X + 2T(p(y)), v(y) = 0. (1.3) 
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Here, p(y) is the inverse ofy{p), which is related to the stream function by p(y) = —^{y), 
while c> is the speed of propagation. The corresponding water depth is given by 

"^Ilvmm- (L4) 

Theorem 2 (Existence of small-amplitude periodic solutions). Let k > be a given wave 
number and assume that the vorticity function is sufficiently regular. Then there is a curve 
of small- amplitude ^--periodic solutions if and only if the Sturm- Liouville problem 

f3M" + - 13' M' = k 2 M, M(p ) = 0, A 3/2 M'(0) = gM(0) (1.5) 

has a nontrivial solution M(p) for some X > — 2r min . Here, (3(p) is defined by 

(3(p) = \ + 2T(p) (1.6) 
for each p < p < 0, while g > denotes the acceleration due to gravity. 

Throughout this paper, we use primes to denote derivatives with respect to p — —ip. The 
main difficulty in finding a necessary and sufficient condition for the solvability of (jl.5p is 
that the solution of this problem is not generally explicit. A typical exception to this rule 
is the case of constant vorticity, in which case (11.51) can be transformed into an ODE with 
constant coefficients. To the best of our knowledge, no other exceptions seem to be known. 

The main goal of this paper is to provide a wide range of non-constant vorticity functions 
for which the eigenvalue problem (jl.5p can be solved explicitly. Our approach is inspired by 
the case of constant vorticity, however, we are able to treat a much larger class of examples 
including polynomial, exponential, trigonometric and hyperbolic vorticity functions. 

Our main result is Theorem [3l which is given in Section [2] together with its proof. Here, 
we use a rather technical change of variables, yet this is inspired by the case of constant 
vorticity that we briefly recall in the beginning of the section. Several applications of our 
main result are given in Section [31 where Examples H] through [TT] are discussed. 

The dispersion relation provided by Theorem [3] appears in (12.1 3p and obviously applies for 
all examples, so we shall not bother to restate it in each case. We only do so in the cases of 
linear vorticity (Example [6]) and exponential vorticity (Example [TT]) . as those are cases for 
which the dispersion relation (I2.13P can be simplified quite a bit. 



2. The main result 



In this section, we present our main result regarding explicit solutions of 

0(p)M"(p) + ~ , (p)M'(p) = k 2 M(p), (2.1) 

where (3(p) = A + 2T(p) is given and k > is the wave number. Our first step is to express 
this equation in terms of the physical variables, see also [21 Lemma 2.3]. We note that our 
definitions ([rT]) - (0]) imply 

^ = v / A + 2r(p(y)) = c - u(y), ^ = -T'(p(y)) = - 7 (-p(y)) (2.2) 

as well as p(d) = and p(0) = po. Using the change of variables 

V (y) = (c-u(y))-M(p(y)), (2.3) 
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we can thus transform equation (12.11) into the Rayleigh equation 

m - *V<») = (2.4) 

u{y) - c 

where dots denote derivatives with respect to y. We now recall (12. 2p and express u{y) in 
terms of the vorticity. This allows us to write the Rayleigh equation (12. 4 p as 

(p{y) - k 2 ip(y) = -i(-p{y))tp(y). (2.5) 

If the vorticity is linear, then this equation is linear with constant coefficients, so it can 
certainly be solved explicitly. Assume the vorticity is constant for simplicity. If we impose 
the boundary condition M(p ) = as in (II. 5p . then (p(0) = by (12. 3 p and so 

<f(y) = c x sinh(%) 

for some constant c\. The corresponding solution of the original problem (12. ip is 

= ci sinh(fcy(p)) = c x ^ / f p kds \ 

{p) ~ v/a+2i» " ^/m' sm [loVm)' 

where P(p) = A + 2T(p), and every solution of (12. ip with M(po) = has this form. 

The computation above shows that the solution of the eigenvalue problem ( 11.5D can be 
determined explicitly, if the vorticity is linear. To extend this observation to more general 
vorticity functions, we now look for solutions of (I1.5P that have the form 



M(p) = J^.sinhH( P ) (2.7) 

for some functions G, H. As it turns out, the unknown functions G, H must be related in a 
very simple way. More precisely, (12. 7p satisfies the ODE in (ll.5p if and only if 

[2PGG" - f3{G') 2 + p'GG' - 2/3"G 2 + A(3G 2 {H') 2 - 4k 2 G 2 } sinhtf 

+ [A(5G 2 H" + A(3GG'H' + 2f3'G 2 H'} cosh H = 0. (2.8) 
To ensure that the coefficient of cosh H is zero, we have to ensure that 

H" G' f3' 
IT + ~G + 2/3 ~ 
and this can be trivially integrated to yield the equivalent condition 



H'(p)G(p)^m = a (2.9) 
for some a e M. Assuming (I2.9p and returning to (I2.8p . we end up with the ODE 

2/3GG" - P{G') 2 + p'GG? - 2f3"G 2 = Ak 2 G 2 - 4a 2 . (2.10) 
As an ODE for 0, this can be easily solved by writing the left hand side as 

2(5GG" - (5{G') 2 + P'GG' - 2p"G 2 = -2G 3/2 ^G 3/2 (^j ^ 

and then integrating twice. This approach already generates a large number of examples, 
but one can also replace the hyperbolic sine in (12. 7p by a regular sine, in which case (12. 9 j) is 
still applicable and the only change in (I2.10p is that —4a 2 becomes 4a 2 . We shall actually 
treat these cases together in our next (and main) result. 
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Theorem 3. Let po < < k. Suppose (3(p), G(p) are positive and C 2 on [po, 0] with 



G 3/2 



2z - 2k 2 G 2 
G^ 2 



for some constant zGl Define the bifurcation parameter A and the vorticity 7 by 

A:=/3(0), 7W=7(-P)==^(P) 



(2.11) 



(2.12) 



and consider the eigenvalue problem ( 11. 5ft which corresponds to this vorticity. Then ( 11. 5ft 
solvable if and only if the dispersion relation 

9 G'(0) 7 (0) 



+ 



A 3 / 2 2G(0) A A 1 /2G(0)F(p ) 
holds, where F(p ) is given by 



(2.13) 



^(Po) = < 



-^tanhf v ~ zd " ifz>0 
f° ^7= i/ z = 



(2.14) 



There are several things that are worth noting about this theorem. 

1. Given any function G, equation ( 12.1 ip can be easily integrated twice to obtain /3, and 
this also determines the vorticity 7 by means of ( I2.12p . Given any function G, in particular, 
one can find a vorticity function 7 for which the dispersion relation is given by (12. 131) . 

2. In practice, one would like to determine the function G in terms of the vorticity 7 
and not the other way around. This is generally a difficult task because it involves solving a 
nonlinear instead of a linear ODE. Nonetheless, it is still feasible for a wide range of vorticity 
functions, as we show in Examples H] through [TTJ 

3. The integral ( I2.14p that appears in the dispersion relation is an elementary integral 
only for two of our eight examples; see the case of linear vorticity (Example [6]) and that of 
exponential vorticity (Example [TTj) . 



4. Our assumption that f3 is positive is equivalent to the assumption A > — 2r min that is 
needed for the trivial flows of Lemma [1] to be defined. Our assumption that G is positive 
can easily be replaced by the assumption that G has one sign; namely, equation ( 12.1 ip is 
equivalent to equation ( I2.18P and G satisfies the latter if and only if — G does. In fact, our 
assumption on G is not method-driven; it is closely related to the fact that the eigenvalue 
problem (II. 5p describes the associated ground state, see pp. 

5. The bifurcation parameter A that appears in the dispersion relation (I2.13P can also be 
expressed in terms of the shear flows of Lemma HJ namely vA = c — u(d) by (ll.3l) - (ll.4p . 

6. The dispersion relation (I2.13P is the necessary and sufficient condition for the solvability 
of the eigenvalue problem (II. 5p . The actual solution of this problem appears in (I2.20p for 
the case z > 0, in (I2.2ip for the case z < and in (12.221) for the case z — 0. 

Proof. We temporarily ignore one of the boundary conditions in (11.51) and focus on 



f3M" + g P' M ' = k2M > M (Po) = 0. 
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(2.15) 



To solve this explicitly, we recall that /3, G are positive and we change variables by 



M{p) = yW)'* ip) - (2 ' 16) 

Then a short computation shows that (12.151) is equivalent to 
4/3GV + 2{2(5GG' + f3'G 2 )^' 

+ (2f3GG" - f3{G') 2 + 0'GG' - 2/3"G 2 - Ak 2 G 2 ) V = (2.17) 
subject to ty(po) = 0. Note that our assumption (12. lip can also be written as 

2/3GG" - P{G') 2 + j3'GG' - 2f3"G 2 - Ak 2 G 2 + Az = 0. (2.18) 
This allows us to reduce the problem (I2.15P to the equivalent problem 

2/3GV + (2(5GG' + f3'G 2 )^' - 2z^ = 0, (p ) = 0. (2.19) 
Suppose first that z > 0. Then the general solution of the last ODE is given by 

ty(p) = C\ sinh I / — — \ + C2 cosh ^ ^ V 



Pa 



and the boundary condition ^f(po) = requires that C2 = 0. In particular, every solution of 
the problem ( 12.15ft has the form 



for some constant C\. We now turn to the original eigenvalue problem (11. 5p . namely (I2.15P 
with an additional boundary condition which reads 

M'(0) _g_ 
M(0) A 3 / 2 ' 

Using logarithmic differentiation, we see that (I2.20p satisfies this condition if and only if 

g G'(0) P'(0) + 1 



A 3 / 2 2G(0) 2/3(0) G(0) v ^(0)F( J 9o)' 

where F(p ) is given by (|2H]) . Since /3(0) = A and /3'(0) = 27(0), the result follows. 
When z < 0, our previous approach applies with minor changes to give 



instead of (I2.20p . while the analogous solution for the case z = is 



W «P) 7 W G(s)i/f)(7j 
In either of these cases, the dispersion relation can be obtained exactly as before. ■ 



3. Some examples 

In this section, we give several applications of Theorem |3j More precisely, we give several 
choices of positive functions /3, G such that (12. lip holds. In each case, the dispersion relation 
is given by (I2.13P and we shall only bother to restate it for Examples M and [TT] namely the 
ones for which the dispersion relation can be simplified to a large extent. 

In our first two examples, we start with a function G and find a vorticity function 7 for 
which our theorem applies. In all other examples, we follow the exact opposite approach 
and start with the vorticity function, instead. 

Example 4 (Polynomial vorticity). Let G(p) = (1 — ap) n for some a > and some real 
number n other than ±2,2/3. Integrating (12. lip twice gives 

m = 2 Cl(1 - apr + 2e2(1 - ap) r /2 + 4:{ !- ap V + 4kHl ~ ap? 



a(3n-2) a 2 (3n-2) 2 a 2 {n 2 - 4) 

The coefficients ci,c%,z can be arbitrary, but we do need /3{p) to be positive on [po,0]. An 
easy way to ensure this is to assume c\,C2,z are non- negative and n > 2. The vorticity 7 is 
then determined using ( 12.121) . In the special case c 2 = z = 0, one finds that 



7(V0 = — onci(l + gm/>) 



n _, 4/c 2 (l + a^) 



a(n 2 -4) ' 

where ^ = — p as usual; the only restrictions in this case are a > 0, c\ > and n > 2. 

Example 5 (Exponential vorticity #1). Let G(p) = e ap for some a ^ and let z e R be 
arbitrary. Then (12.111) can be easily integrated twice to give 

Ab 2 A 7 p~ 2a P 

Pip) = + c 2 e-^ 2 + % + ^-5-. 

If we require ci,c%,z to be non- negative, then /3 is positive and our theorem applies. Once 
again, the corresponding vorticity is given by (12. 12ft . namely 

ac x ac 2 e a ^l 2 Aze 2a ^ 
Example 6 (Linear vorticity). Suppose 7(V0 = aip + b for some a, 6 G M. Then 



Pip) = A 



/■p 

+ 2 / (b - as) ds = X + 2bp - ap 2 
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so 0"{p) = —2a. In particular, (12. lip holds with G(p) = 1 and z = k 2 — a. The dispersion 
relation (12. 13[) provided by Theorem [3] involves the integral 

ds f° ds 



Po G(s)y^{s) J P0 ^X + 2T(s) 

which is equal to the depth d by (II ,4p . Thus, the dispersion relation ( 12 . 1 3[) reduces to 

A - bF(d)X 1/2 - gF(d) = 0, (3.1) 
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where F(d) is given by (I2.14p . namely 

{tanh(d\/ k 2 — a) 
Vk 2 -a 
d 
tan(d\/<i— k 2 ) 
Va-k 2 

Needless to say, one can also solve the quadratic equation (13.11) and write it as 

>/A = ^ ± \ y/PF(d)* + AgF(d), (3.2) 

where the sign is chosen so that v^A = c — u(d) is positive. To the best of our knowledge, 
this dispersion relation was only known in the case of constant vorticity, in which case it 
goes back to [3] , see also [TJ [3] • 

Example 7 (Quadratic vorticity). Suppose j(if}) = aip 2 for some a^O and let 

G(p) = - 20a 2 k 2 p 2 + lQak^p + (25a 2 A - 16fc 6 ), f3(p) = A + 

Then one can verify that (12.181) . or equivalently (12.111) . holds with 

z = A; 2 (1125a 4 A 2 - 1056a 2 fc 6 A + 256k 12 ) > 0. 

We actually came up with this example by looking for solutions of (I2.18P for which /3, G are 
both polynomials. The same approach applies to any quadratic vorticity whatsoever, but 
the corresponding cubic solution G looks much messier in the general case. Note that is 
monotonic for this example and that the same is true for G, since 

G'(p) = 2a(25aV - 20ak 2 p + 8k A ) 

and the discriminant of the quadratic is — 400a 2 fc 4 . In particular, positivity of 0, G on the 
interval \po, 0] amounts to positivity of /3, G at the endpoints of the interval. 

Example 8 (Cubic vorticity). Suppose "f(ip) = atp 3 + hip for some a, b G R and let 

G(p) = + 3a(k 2 + 3b)p 2 + (2k 4 + 6bk 2 - 9aA), (3(p) = A - — - bp 2 . 

Then one can verify that (12.181) holds with 

z = 2k 2 {2k 4 + 6bk 2 - 9a\)(k 4 + 2bk 2 - 6a\ - 3b 2 ). 

This example was also obtained by looking for polynomial solutions of f !2. 18j) . an approach 
that no longer works for a more general cubic or even a quartic vorticity. To ensure that our 
theorem is applicable on any interval [po, 0] whatsoever, we impose the conditions 

G(0) = 2k 2 (k 2 + 3b) -9a\ > 0; a < < A; b < 0. (3.3) 

Assuming those, (3 is decreasing with (3(0) > 0, so it is positive on [p , 0], while 

G'(p) = 6ap(3ap 2 + k 2 + 3b). 

If it happens that k 2 + 3b < 0, then G is also decreasing with G(0) > 0, so we are done. If 
it happens that k 2 + 3b > 0, then G has a unique critical point p* < 0, where 



Jk 2 + 3b 3 
G{p.) = -9a\ + Uk 2 + 3b)(k 2 -b)>0. 
3\a\ 2 



In either case then, G is positive on the whole negative real £1X1 S ; cLS needed. To avoid the 
technical condition in (I3.3p . one may focus on the special case a < < A and 6 = 0. 

Example 9 (Trigonometric vorticity). Suppose ^(ip) = aicos(bip) + c^sin^^) for some 
ax, a,2 6 M and some b > 0. In this case, one can take 

G( P ) = m + ^f, m = a + 2aiSi ^ bp) + 2 <<*»m - ^ 



b 2 ' ^ y b b 

and check that (I2.18p holds with 

k 2 



b 4 



(Ak 2 + \b 2 ) 2 - Aa 2 b(Ak 2 + Xb 2 ) - Aa\b 2 



To ensure that (3(p) is positive for all p, we have to impose the condition 



2 



A > - ( a 2 + yal + at 

while the positivity of G(p) follows trivially from that of f3(p). This example is related to 
the previous two in the following sense. When a\ = 0, the eigenvalue problem (11. 5p can be 
solved using the substitution u = cos(6p) that reduces the problem to one with polynomial 
coefficients; the solution G(p) was found by looking for polynomial solutions in u. 

Example 10 (Hyperbolic vorticity). Suppose ^(ip) = a± cosh(bip) + a 2 sinh(6-?/') for some 
ax, a,2 6 M and some b > 0. In this case, one can take 

G{p ) = fKp)-% m = x+ 2a 1 sinHb P ) + 2a 2 (l-cosHb P )) 
b z b b 

and check that (I2.18P holds with 

z = ^r- \(4k 2 - Xb 2 ) 2 - Aa 2 UAk 2 - Xb 2 ) + Aa 2 b 2 } . (3.5) 
o 4 L J 

Needless to say, one can always obtain this example from the previous one by allowing the 

coefficients to be complex. If we assume ai < min{— a 2 , 0}, then we have 

(5\p) = (ai - a 2 )e bp + (at + a 2 )e~ bp < 

for all p < 0. This is clear when a\ — a 2 < 0, in which case (3' is the sum of two non-positive 
terms, but it also holds when a\ — a 2 > 0, in which case 

0'(p) < (ai - a 2 )e~ bp + (a, + a 2 ) e - bp = 1a x e~ bp < 0. 

It easily follows that our theorem applies when a\ < min{— a 2 , 0}, b > and A > 

Example 11 (Exponential vorticity #2). Suppose ^(ip) = ae b ^ for some a, b ^ 0. This is a 
special case of the previous example and one can take 

, , . Ak 2 . . , 2a(e- bp -l) k 2 (Ak 2 - Xb 2 - lab) 2 . , 
G(p) = /3(p) - — , /3(p) = A-— ^ z = ^ L (3.6) 

In this case, the dispersion relation (I2.13P can be simplified quite a bit. Suppose z > for 
the moment. Then (I2.20p shows that every solution of (11.51) is a scalar multiple of 



Mo(p) = W f^' Sinh 




hence also a scalar multiple of 



. r G(p) . , / f p kUk 2 - Xb 2 - 2ab) 7 . 



We now recall our choice (I3.6P of G(p) and we make use of the identity 

k k{Ak 2 - Xb 2 - 2ab) 



i, ^W) + zk/b 

- log - !L 1 = 

2 y/ftp)-2k/b 



(3.8) 



that may be viewed as a partial fractions decomposition for the right hand side. Using this 
identity and our choice of G(p) in (I3.6p . one finds that 



M M = ( i+ -£= ) J^R- 2k ' b . ex p ' r 



i ,, ^ IVPGS + n/i exp r _^ 




, 2 t>VWF)J y VPM-2k/b y„vW, 

Define the coefficients Ci, C2 by the formula 

ci = v 7 ^) - 2fc/6, c 2 = y^) + 2A;/6. 
In view of (II. 3p . these can also be expressed in terms of the velocity field as 

Cl = c-u(0) -2k/b, c 2 = c-u{0) + 2k/b. (3.10) 
Now, every solution of (11.51) is a scalar multiple of (13.91) . hence also a scalar multiple of 

M(p) = ci 

I o u \ I rv u a „ \ 

, , ■ - - 1 ' . , , , . (3.11) 

We remark that the last formula was obtained under the assumption that z,/3,G are all 
positive, yet this formula provides a solution of (II. 5p as long as merely (3 is positive. Next, 
we turn to the boundary condition at p — 0, which reads 

A 3/2 M'(0) = gM(0). (3.12) 

Differentiating ( 13. lip , one finds that 

kci A 2k (3'(p)\ 
M'{v) = =L= 1 + — == + f-f^ • exp 

kc 2 ( 2k ?(p)\ 

b^m WW) ^ 
In view of (13. 6p and ( II. 4p . this implies 

kc x e~ kd ( 2k 2a\ kc 2 e kd ( 2k 2a 

M (0) = 1 - 7 =^- 1 + — = + — V- 1 = + — 

VX \ bVX bXJ y/\ V bVX bX 




and also 



M(0) = ae 



-kd 



1 + 



2k 



c 2 e 



kd 



1 - 



2k \ 



bsfXJ 



■ 



Thus, the boundary condition (I3.12p holds if and only if 



-kd 



(g + kX) 



fi 



2k 



+ 



2ak 



c 2 e 



kd 



(g -kX)(l 



2k 



2ak 



Multiplying through by byX, one thus obtains the dispersion relation 

c x ({g + k\){bV\ + 2k) + 2ak^ = c 2 e 2kd ((g - k\)(bVX — 2k) - 2akV\), (3.13) 
where d is the depth, Va = c — u(d) and the coefficients ci, c 2 are given by f)3.10p . 
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